Astronomy & Astrophysics manuscript no. 16899 © ESO 201 1 

July 5, 2011 



^T 



o 



Multi wavelength campaign on Mrk 509 

III. The 600 ks RGS spectrum: unravelling the inner region of an AGN 

R.G. Detmers 1 ' 2 , J.S. Kaastra 1 - 2 , K. Steenbrugge 3 ' 4 , J. Ebrero 1 , G. A. Kriss 5 ' 6 , N. Arav 7 , E. Behar 8 , E. Costantini 1 , G. 
Branduardi-Raymont 9 , M. Mehdipour 9 , S. Bianchi 10 , M. Cappi 11 , P. Petrucci 12 , G. Ponti 13 , C. Pinto 1 , E.M. Ratti 1 , and 

T. Holczer 8 

1 SRON Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands e-mail: 
r . g . detmers@sron . nl 

2 Astronomical Institute, University of Utrecht, Postbus 80000, 3508 TA Utrecht, The Netherlands 

3 Instituto de Astronomfa, Universidad Catolica del Norte, Avenida Angamos 0610, Casilla 1280, Antofagasta, Chile 
f"^ ' 4 Department of Physics, University of Oxford, Keble Road, Oxford OX1 3RH, UK 

{SJ ! 5 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

Department of Physics & Astronomy, The Johns Hopkins University, Baltimore, MD, 21218, USA 
' Department of Physics, Virginia Tech, Blacksburg, VA 24061, USA 

Department of Physics, Technion, Haifa 32000, Israel 

Mullard Space Science Laboratory, University College London, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK 

Dipartimento di Fisica, Universita degli Studi Roma Tre, via della Vasca Navale 84, 1-00146 Roma, Italy 

INAF-IASF Bologna, Via Gobetti 101, 1-40129 Bologna, Italy 

UJF-Grenoble 1 / CNRS-INSU, Institut de Plantologie et d?Astrophysique de Grenoble (IPAG) UMR 5274, Grenoble, F-38041, 

France. 
\,J ■ 13 School of Physics and Astronomy, University of Southampton, Highfield, Southampton S017 1BJ, UK 

C^; abstract 

O ' 

5_l ■ We present the results of our 600 ks RGS observation as part of the multiwavelength campaign on Mrk 509. The very high quality of 

the spectrum allows us to investigate the ionized outflow with an unprecedented accuracy due to the long exposure and the use of the 

RGS multipointing mode. We detect multiple absorption lines from the interstellar medium and from the ionized absorber in Mrk 509. 

A number of emission components are also detected, including broad emission lines consistent with an origin in the broad line region, 

the narrow O vn forbidden emission line and also (narrow) radiative recombination continua. The ionized absorber consists of two 

velocity components (v = - 13 ± 1 1 km s -1 and v = — 319 ± 14 km s~'), which both are consistent with earlier results, including UV 

data. There is another tentative component outflowing at high velocity, - 770 ±109 km s" 1 , which is only seen in a few highly ionized 

absorption lines. The outflow shows discrete ionization components, spanning four orders of magnitude in ionization parameter. Due 

l/~) ■ to the excellent statistics of our spectrum, we demonstrate for the first time that the outflow in Mrk 509 in the important range of log £ 

\& • between 1-3 cannot be described by a smooth, continuous absorption measure distribution, but instead shows two strong, discrete 

^^ ' peaks. At the highest and lowest ionization parameters we cannot differentiate smooth and discrete components. 
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• • , 1. Introduction 2000). While we have a reasonable qualitative understanding of 

*•*■ ■ the feedback from relativistic jets ( as observed in clusters of 

V-| ■ One of the main reasons to study active galactic nuclei (AGN) galaxies , see e.g. Fabia nlTalll^OOl) . we still lack a quantita- 

is to learn about feedback from the AGN to the galaxy and its tive picture of the fee dback of the AGN on the galaxy and on its 

S . direct environments. Feedback is a combination of enrichment surroundings 

' (the spreading of elements into the interstellar and inter galactic There is a broad ongoing effort to improve this, and recent work 
media (ISM and IGM), momentum feedback (due to winds), and on broad absorption i ine (BAL) quasars shows that the mass out- 
direct kinetic feedback (i.e. energy ejection into the ISM or IGM flow rates in these systems w& 10 0s of solar masses per year 
by jets). From recent observations on cooling clusters of galax- and the kinetic i urmnosity involved is a fe w percent of the to- 
les (see e.g. McNa mara & Nulserj | 2007j for an overview), as tal bolometric luminosity dMoe et alj r2009: Pu nnet alJl20Tf)h . 
well as from recent insights into galaxy and AGN co-evolution There have also been i ndications that some A GN harbor a highly 
dDiMatteoeta l. 2005; Elvis 2006; Bower 2009; Fab iarj [2010t), i onize d, massive, ultra-fast o utflow, wi th velocities reaching up 
it has become clear that feedback from AGN is a crucial ingre- to 60000 km s' 1 dReeves et alJ|2003b iPounds & Ree ves 2009; 



dient the evolution of galaxies and clusters of galaxies. This is |p ont j e t all l2009t iTombesiet alJEnOaTlbl) . These outflows are 

also seen in the so-called M-cr relation, which links the velocity hard t0 detect , howeV er, and appear t0 be var i a ble (because they 

dispersion of stars in the bulge to the mass of the Super-Massive &K only present in some obserV ations of a single source). These 

Black Hole (SMBH) flFerrarese & Merntt||2000|; |Gebhardt et alj m extreme cases of outflows that are present in only a frac- 

tion of the total number of AGN. Whether feedback from less 
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extreme outflows, such as those that are present in about 50 % 
of the local Seyfert 1 galaxies is also important, remains an un- 
solved question. 

If we can establish the impact that the outflow has on the 
galaxy in these local Seyfert 1 galaxies, this should allow us 
to extend the feedback estimates that we obtain to higher red- 
shifts to the more powerful AGN, which we are unable to in- 
vestigate with the current generation of X-ray grating spec- 
trometers. However we first need to deal with the two main 
uncertainties concerning the outflows. These are the geometry 
of the inner region o f an AGN and the locat i on or origin of 



the o utflow (see e.g 
l2001t IdaskelletalJ 



iMurrav & Chiang] 119971; JKrolik & Krissl 
2007t) . Earlier work has placed the out- 
flow at various distances and also the est i mates for feedback 
can vary wildly (see e.g. JBehar et al]|2003t iBJustin et alj|2005t 
Krongold et al. 2007; Detmers et al. 2008, for some examples). 
Therefore answering these two questions is the main goal of the 
Mrk 509 multiwavelength campaign. 

Multiwavelength campaigns on AGN are crucial for gaining 
a complete understanding of the inner regions of these sources. 
Earlier multiwavelength campaigns focus ed mainly on a bun- 
dance studies of the outflow (see Mrk 279, Arav et aLll2007l) or 
on determining the outflow structure and location by combin- 
ing UV and X-ray data eg. NG C 5548 ( Steenbrugge et al. 2005); 
NGC 3783 (iNetzer et al.l2003LlGabel et al.l20 03). Our dedicated 
multiwavelength campaign on Mrk 509 is much more exten- 
sive than previous attempts. Our more intensive observations are 
ideal for locating the outflow, using the variability of the source 
and response of the ionized gas to determine its location (the use 
of variability to locate gas has been very successfully used in 
reverberation mapping of the BLR, see e.g. iPeterson & Wandel 
2000: lDennev et aDl2010l, for an overview of the method and for 
the latest results). 

Apart from the location and kinematics, one of the other im- 
portant questions regarding the outflow is what the ionization 
structure is. Earlier studies have reported different results. The 
outflow in NGC 5548 appear s to be a continuous distribution of 
column density vs. log £, (ISteenbrugge et alJl2005h . NGC 3783, 
on the other hand, shows differe nt separate ionizatio n compo- 
nents, all in pressure equilibrium (Krongold et al. 2003). In Mrk 
279 the situation appears to be more complex, because a non- 
monotonous, contin uous distribution p rovides th e best descrip- 
tion to the data (Costant ini et al.l l2007). Recently. iHolczer et alJ 
(120071) and Behar (2009) have shown that for most local Seyfert 
1 galaxies with an outflow, a continuous distribution of column 
density vs. ionization parameter is the best description of the 
data. What they also show is that there are distinct £ values where 
no is gas present. They interpret these gaps as thermal instabil- 
ities that cause the gas to rapidly cool or heat and then shift to 
other ionization parameters. What is clear from these studies is 
that there is no single model that describes all the observed out- 
flows. High-quality, high -resolution spectra of the outflows are 
crucial for investigating the structure, since it can be the case (as 
in Mrk 279) that some components of the outflow have very low 
column density, which would otherwise escape detection. 

Mrk 509 is one of the best studied local AGN, and due to its 
large luminosity (L(l - 1000 Ryd) = 3.2 x 10 38 W), it is also con- 
sidered to be one of the closest QSO/Seyfert 1 hybrids. Earlier 
work on the outflow in the X-ray regime has revealed that it con- 
sists of a wide range of ionization components, but lacks the very 
high and also very low ionized gas (weak Fe UTA and no Si xiv 
Lva. I Yaqoob et al. 2003). The outflow has been described using 
three i onization compon ents, each with a different outflow ve- 
locity (ISmifh et al.l l2007); however, the exact outflow velocities 



differ between differe nt publications, most l ikely due to a limited 
signal-to-noise ratio. lDetmers et alJ (12010) have analyzed three 
archival observations of Mrk 509 with XMM-Newton. They also 
find three com ponents for the outflow, although with different 
velocities than Smith et al. (2007). Including the EPIC-pn data 
and improving the relative calibration between RGS and EPIC- 
pn achieved increased sensitivity. With this improvement they 
were able to detect variability in the highest ionization compo- 
nent, constraining the distance of that component to within 0.5 
pc of the central source. Another point of interest is that there 
have been indications o f an ultra-high velocity outflow a s seen 
through the Fe K line (ICappi et al.ll2009HPonti et al.l|2009l) . This 
outflow could make a potentially large contribution to feedback, 
as the velocity is very high, although it appears to be transient 
dPonti et al.ll2009h . 

This work is the third in a series of papers analyzing the 
very deep and broad multiwavelength campaign on Mrk 509. 
The co mplete campaign details are presented in Kaastr aet alJ 
(201 la), hereafter Paper I. Here we present the main results ob- 
tained from the stacked 600 ks XMM-Newton RGS spectrum. 
With this spectrum, we are able to characterize the properties of 
the ionized outflow in great detail (velocities, ionization states, 
column densities, density profile, etc.). Other features detected 
in the spectrum ( emission lines, Galactic absorption, etc.) will 
not be discussed in detail here. Because different physics are in- 
volved, we will discuss them in future papers in this series. 

This article is organized as follows. Section [2] briefly de- 
scribes the data reduction for obtaining the stacked spectrum, 
and we show the spectral models that we use to describe the 
data. The spectral analysis and the results are presented in Sect. 
|3] We discuss our results in Sect.|4]and present our conclusions 
in Sect. |5] 



2. Data reduction and modeling 

2.1. Data reduction 

The RGS data reduction used here is much more complex 
than the standard pipeline processing using the XMM-Newton 
Science Analysis System (SAS), the main reasons among others 
being the use of the RGS multi-pointing mode, a variable source 
and a nonstandard procedure of filtering for bad pixels. The full 
d etails of the data reduc tion and all the necessary steps are found 
in Kaastra et alJ (l2011bl) . hereafter paper II. In short, we used the 
SAS 9.0 software package to reduce all the individual observa- 
tions. Then we created a fluxed spectrum for each observation 
and stacked those taking the effects of the multi-pointing mode 
into account. This way a fluxed RGS spectrum was created by 
stacking both RGS 1 and 2 and both spectral orders. We used 
this fluxed spectrum for fitting our data. Figure Q] shows the full 
fluxed spectrum with some of the strongest lines indicated. 

2.2. Setup 

We used data between 7 and 38 A when fitting the RGS spec- 
trum. As the fluxed spectrum consists of both the RGS 1 and 2 
data in both spectral orders, we binned the data between 7 and 38 
A in 0.01 A bins. The average number of counts per 0.01 A bin 
is approximately 900 (i.e. a signal-to-noise-ratio of about 30). 
We therefore used^ 2 statistics when fitting the spectrum. All er- 
rors are given for A^ 2 = 1. We used the SPEX 2.03.00J spectral 
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Table 1. Spline continuum parameters. 



20 
Wavelength (A) 

Fig. 1. The fluxed stacked RGS spectrum in the 7 - 38 A range. 
The strongest lines are indicated and the Oi ISM edge can be 
clearly seen around 23 A. 



fitting package to fit the spectrum. We updated the wavelengths 
of some important transitions for our study (see AppendixlAl. 

We constructed the spectral energy distribution (SED) for 
Mrk 509, using the EPIC-pn and OM data to obtain the neces- 
sary flux points for the XMM-Newton observations and extend- 
ing it with other data. This SED is an average of the Mrk 509 
SED during the time of observations (roughly two months time). 
Th e full procedure on ho w the SED was derived can be found 
in iKaastra et alJ (l2011al) . The ionization balance calculations 
needed for our spectral modeling (the xabs components, see 
Sect. 13. 2t were bas ed on this SEP and p erform ed using version 
C08.00 of Cloud\PI (lFerland et alJll998h with lLodders & Palmd 
(2009) abundances. 

2.3. Spectral models 

The unprecedented quality of the spectrum requires a rather 
complex spectral model to be described accurately. We describe 
each component in more detail in separate sections, but we give 
a short overview of the total model here. 

We model the continuum with a spline (see Fig.|2|. The main 
reason for doing so is that a spline can accurately describe the 
(complex) continuum shape without having to make any physi- 
cal assumptions about the origin of the shape of the continuum 
(powerlaw, blackbody, Comptonization, or reflection, for exam- 
ple). We use a r edshift z - 0.03450, which combines the cosmo- 
logical redshift dHuchra et al.ll 19931) with the orbital velocity of 
the Earth around the Sun, which is not corrected for in the stan- 
dard XMM-Newton analysis (see Kaastra et al. 201 la). Galactic 
absorption (N H = 4.44 x 10 24 irT 2 , Murph y etalJI 19961) is also 
taken into account. We use three distinct phases for the Galactic 
ISM absorption, a neutral [kT - 0.5 eV) phase, a warm (kT = 4.5 
eV) slightly ionized phase, and a hot (kT - 140 eV) highly ion- 
ized phase (tPinto et al.l2010l) . Additionally we model the neutral 
oxygen and iron edges of the ISM by including a dusty compo- 
nent. Details about the Galactic foreground absorption are given 
by Pinto et al. (in prep). 

The ionized outflow is modeled with three different mod- 
els, each with multiple (two or three) velocity components to 



Wavelength 


Flux 


(A) 


(photons m 2 s _1 A -1 ) 


5.00 







5.95 


0.7 


±0.6 


7.07 


13.95 


±0.09 


8.41 


14.28 


±0.06 


10.00 


15.10 


±0.04 


11.89 


15.92 


±0.09 


14.14 


17.27 


±0.09 


16.82 


19.36 


±0.16 


20.00 


21.31 


±0.11 


23.78 


25.15 


±0.08 


28.28 


27.83 


±0.11 


33.64 


33.62 


±0.19 


40.00 


9.80 


±0.18 



account for the separate outflow velocities observed. All mod- 
els take a wide range of ionization states into account. These 
models are described in more detail in Sect. 13.21 We also in- 
cluded eleven broad and narrow emission lines, which are mod- 
eled with a Gaussian line profile. Radiative recombination con- 
tinua (RRCs) are also included using an ad-hoc model that takes 
the characteristic shape of these features into account . 



3. Spectral analysis 

3.1. Continuum, local absorption, and emission features 

The continuum is modeled with a spline with a logarithmic spac- 
ing of 0.075 between 5 and 40 A. We show the spline model in 
Fig. [2] and in Table Q] The continuum is smooth, so the spline 
does not mimic any broad line emission features. The softening 
of the spectrum at longer wavelengths can be seen clearly. 
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20 
Wavelength (A) 

Fig. 2. The unabsorbed spline continuum model used for the 
Mrk 509 observations. 



The neutral Galactic absorption is responsible for the narrow 
O i (23.5 A) and N i (3 1 .3 A) absorption lines. To fit the Galactic 
O ii absorption line we add a slightly ionized component with a 
temperature of 4.5 eV and with a column density that is 4 % of 
the cold, neutral gas (Pinto et al., in prep). To properly model the 
O i edge, we take the effects of depletion into dust into account. 
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Table 2. Broad emission line parameters for the combined spec- 
trum, where fluxes are corrected for Galactic and intrinsic ab- 
sorption. 



Table 3. RRC parameters. 



Line 



Wavelength 

(A) 



Flux 
(ph irT 2 s _1 ) 



4r 2 



Ovn r 


21.602 


1.00+0.09 


82 


Ovn i 


21.802 


<0.56 





Ovnf 


22.101 


<0.7 





O vm Lyo- 


18.967 


0.42±0.04 


46 


N vii Lytt 


24.780 


0.09±0.06 


1 


C vi Lyo- 


33.736 


0.25±0.10 


3 


Neix r 


13.447 


0.26±0.03 


40 



All the wavelengths were fixed to the laboratory wavelength. 



This same depletion is also responsible for the lack of a strong 
neutral iron absorption blend around 17.5 A. We use dust that 
consists of pyroxene and hematite. 

The sightline to Mrk 509 passes thro ugh a high-velocit y 
cloud located in the Galactic halo (see e.g. lSembach et alll 995). 
The hot gas as seen in the C vi, O vn, O vm, and Ne ix absorption 
lines at z = could either be from ISM absorption in our Milky 
Way or from this high-velocity cloud (which has an LSR veloc- 
ity of around - 250 km s~'). A more detailed multiwavelength 
analysis and comparison of all these local components will be 
presented in a follow-up paper (Pinto et al. ,in prep). 

The spectrum shows some emission lines, most of which are 
broadened (Table |2). In this table we also show the change in 
X 1 when the line is omitted from the model. The broad emission 
lines are visible as excesses on both sides of the corresponding 
absorption lines (see Fig. [3). In our modeling we fix the width 
of the lines to an FWHM of 4200 + 200 km s _1 as measured 
for the Balmer lines si multaneously by the OM optical grism 
(Medhipouret al. 2011). We assume that these lines originate in 
the BLR (Costantin retal.ll2007l) . 

We also detect the narrow Ovn f emission line with an 
absorption-corrected flux of 0.46 + 0.06 ph m~ 2 s , as well 
as a weaker intercombination line with a flux of 0.13 + 0.04 
ph irT 



■ 2 s ' and an Ne ix f emission line with a flux of 0.09 



0.02 ph m s . The narrow resonance line cannot be sepa- 
rated from the absorption line, so we have assumed a ratio of 
3:1 for photoionized gas for the forbidden-to-recombination 
line ratio (Porquet & Dubau 2000). We do not detect any other 
narrow emission lines. In our spectral modeling we adopt an 
FHWM of 610 km s -1 for these narrow emissio n lines, based 
on th e width of the optical [Oin] emission line (Phillip s et al.l 
1983). We assume here that these lines are produced in the NLR 
(Guainazzi & Bianchi 2007). Radiative recombination continua 
(RRC) have been detected in other Seyfert 1 spectra (see 
e.g.N GC 3783 or Mrk 279, iKasoi et al.ll200l ICostantini et all 
2007), so we checked whether they are present in Mrk 509. The 
RRCs are weak and hard to detect (Table [3). We modeled them 
using the rrc model of SPEX and obtain a temperature of 5.4 + 
2.5 eV for the RRC. This low temperature is a direct indication 
of photoionized gas. 



3.2. Ionized outflow 

We use three different models for characterizing the ionized out- 
flow. We describe each model separately and then compare the 
results obtained by all three models. By comparing the results for 
the different models for the ionized outflow, we can investigate 
its ionization structure and density profile. 



Ion 


Wavelengi 


th Flux 




(A) 


ph irT 2 s" 1 


Ovn 


16.78 


0.08 ± 0.04 


O vm 


14.23 


<0.01 


Cv 


31.63 


0.18 ±0.14 


Cvi 


25.30 


0.16 ±0.08 


Nvi 


22.46 


0.41 ±0.18 


Nvn 


18.59 


<0.08 


Neix 


10.37 


0.06 ± 0.03 


Nex 


9.10 


<0.01 



The first model (Model 1 from here on, TableHJl contains two 
slab components. The slab model of SPEX calculates the trans- 
mission of a slab of material with arbitrary ionic column densi- 
ties, outflow velocity v, and r.m.s. velocity broadening <j as free 
parameters. We assume a covering factor of unity for both com- 
ponents. The slab components have a different outflow velocity 
and velocity broadening. However, we assign each ion to one of 
the components based on the ionization parameter £, which is 
defined in the following way: 



€ 



(i) 



where L is the 1 - 1000 Rydberg luminosity, n the hydrogen 
number density of the gas and, r the distance from the ionizing 
source. For Mrk 509 we obtain an ionizing luminosity of L — 
3.2 x 10 38 W from the SED. All ions with log f < 2.1 are as- 
signed to the first component, while the others are in the second 
component. The division is based on a preliminary fit where the 
outflow velocity of individual ions was a free parameter. The first 
component has an outflow velocity v of - 57 + 8 km s -1 and a 
velocity broadening <x = 158 + 5 km s -1 , while the second com- 
ponent has an outflow velocity y = — 254 + 40 km s and a = 
133 ± 30 km s . The ionization parameters are given in Table|4] 
and are those for which the ion fraction peaks for that ion. The 
fit gives x 1 - 3643 for 3109 degrees of freedom (dof). We also 
show in Table [4] the best-fit velocities for individual ions when 
we leave the velocity free compared to the velocities of all the 
other ions. 

The second model (Model 2 from here on, Table |4]i is an 
extension of Model 1 . Instead of separating the ions according 
to their ionization parameter, we include all ions for both ve- 
locity components. Additionally we add a third velocity compo- 
nent to account for a high- velocity (y = - 770 km s _1 ) com- 
ponent (tentatively detected in the Chandra HETGS spectrum, 
see lYaqoob et al.ll2003l) mainly to get the appropriate line cen- 
troid for the Fexxi and Mgxi absorption lines. This model is 
more realistic than Model 1 because it assumes a multiveloc- 
ity structure for every absorption line, which is consistent with 
what has been observed in earlier UV observations of Mrk 509 
dKriss et aLllJOOOt iKraemer et al]|2003l) . The fit gives x 2 = 3589 
for 3070 dof. 

The third model (Model 3 from here on, Tabled is based on 
multiple photoionization components (xabs). In each xabs com- 
ponent the ionic column densities are related through the ion- 
ization parameter £. Free parameters are the hydrogen column 
density Nh, ionization parameter £, r.m.s. velocity width cr, and 
outflow velocity y. Based on the results of Model 2, we start with 
one xabs component for each of the two main velocity compo- 
nents detected. We add extra xabs components until the fit no 
longer improves. A fit with only one xabs component for each 
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Fig. 3. Broad and narrow X-ray emission lines of O vm Ly a (left) and O vn (right). The model without any broad lines is shown 
as the dashed-dotted line, while the model (model 2) with the lines and RRC included (marked as XBLR and n) is shown as a solid 
line. 



velocity results in a fit of^- 2 = 4261 for 3145 dof. Adding an ad- 
ditional xabs component for each outflow velocity improves the 
fit by Ax 1 - 264/4 dof. We use the same outflow velocity and 
r.m.s. velocity for the xabs components of each velocity compo- 
nent. If we add a third pair of xabs components we again im- 
prove the fit significantly resulting in a further improvement of 
Ax 2 - 138/4 dof. As a last step, we leave the outflow velocity 
and the r.m.s. velocity width for each component free. This then 
improves the fit by Ax 2 - 32/4 dof The best fit we obtain this 
way has x 2 - 3827 for 3157 dof. The results for the final fit are 
shown in Table [5] however the component with the highest ion- 
ization parameter (component E2 in Table |5]l shifts to a much 
higher outflow velocity, namely 492 km s . We therefore add 
it to the fast velocity group, so that the slow outflow can now 
be described properly with two xabs components, while the fast 
outflow is described by four components. We label each compo- 
nent according to its ionization parameter (A to E for increasing 
£) with index 1 or 2 for low and high outflow velocity, respec- 
tively. 

As a test we also fit a fourth model, the so-called warm 
model, essentially a power-law distribution of xabs components . 
This is similar to the model used by Steenbrugge et al. (2005) 
in NGC 5548 (model D in that paper) and akin to the absorp- 
tion m easure distribution (AMD) method used bv lHolczer et al.l 
(2007). We first defined a range of ionization parameters, be- 
tween which we fit our model. In our case we fit the model be- 
tween log £ = — 2 and 4. We used a grid of 19 points in order to 
accurately describe the total AMD. At every grid point a value /; 
was determined, which is defined as /; = dNu /dlog^. This way 
we obtained the distribution of Nh versus log £. We did this for 
both velocity components (ignoring the very fast -770 km s 
outflow at the moment). The best fit we obtain has^ 2 = 3822 for 
3120 dof. However due to the correlation between the two warm 
component parameters (especially the factors /;), calculating the 
exact error on every parameter is very difficult. 

We therefore use the warm model only to check whether we 
have missed any ions in our slab fit. With a continuous model 
like warm, weaker lines that otherwise might be hard to detect 
are taken into account. In this way we have identified several 
ions, which are predicted to be present based on the warm model, 
but were not included in the slab fit since they produce only weak 



lines. These are Ne vn and Si x - Si xn. All these ions, however, 
have very weak lines in the Mrk 509 spectrum, and the fitted 
ionic column densities (using the slab model) only yield upper 
limits. We therefore are confident that Models 2 and 3 are accu- 
rate representations of the data. 



Table 5. Parameters for model 3. 



Comp log £* 



Bl 

CI 



0.81 : 
2.03: 



0.07 
0.02 



0.8: 
2.6: 






logU" 



0.1 

0.2 



124: 
193: 



20 
14 



25 ± 
-43 



30 
b20 



-0.73 
0.49 



A2 
C2 
D2 
E2 



-0.14 ±0.13 
2.20 ± 0.02 
2.62 ± 0.08 
3.26 ± 0.06 



0.4 ±0.1 

4.4 ± 0.5 
1.8 ±0.5 
6.3 ± 1.2 



79 d 

29: 

34 h 
37 h 



26 

:6 

19 
19 



- 180 ±41 
-267 ±31 

- 254 ± 35 

- 492 ± 45 



-1.68 
0.66 
1.08 

1.72 



" Ionization parameter in 10 W m. 

* Column density in units of 10 24 m" 2 . 

c r.m.s. velocity broadening in km s 

d Outflow velocity in km s _1 , a negative velocity corresponds to a blueshift. 

' ionization parameter (as used in UV spectroscopy). 



3.3. Spectral fit 

The RGS spectrum and the best-fit model (Model 2) are shown 
in Figs. |4] - |8] All the strongest absorption lines are labeled. 
Galactic lines are indicated with z — 0. The spectrum has been 
re-binned to 0.02 A bins for clarity. The model reproduces the 
data very well. We detect the O vm Lyman series and the O vn 
resonance transitions up to the ls-5p transition, as well as the 
C vi Lyman series up to the ls-6p transition. We also detect the 
O vn and O vm series from the local z = component. The only 
features that are not reproduced well are the N vi ISM absorption 
line at 28.78 A, a feature around 32.5 A, and another near the 
Cv absorption line at 33.9 A. The 32.5 A feature is most likely 
due to small residuals in the RGS calibration, because it is much 
sharper and narrower than the other emission / absorption fea- 
tures. Its wavelength does not correspond to known major tran- 
sitions. Also the O vi absorption line at 22.8 A is not well-fitted, 
possibly due to blending with the O vn f narrow emission line. 
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Table 4. Predicted and measured Ionic column densities for the most important ions of the outflow for model 1 and model 2. 



logf 



Nf 



N" 



:.n:v 



"A? A?" 



D' E^ 



N A 



N 



N'" 



-8.50 

-8.50 

-8.50 

-2.10 

-1.65 

-1.48 

-1.45 

-1.42 

-0.95 

-0.65 

-0.60 

-0.55 

-0.55 

-0.05 

0.00 

0.05 

0.15 

0.42 

0.75 

0.85 

0.92 

1.15 

1.20 

1.22 

1.35 

1.37 

1.60 

1.65 

1.71 

1.75 

1.85 

1.88 

1.97 

2.00 

2.04 

2.05 

2.10 

2.15 

2.20 

2.21 

2.25 

2.30 

2.35 

2.42 

2.42 

2.51 

2.60 

2.77 

3.01 

3.20 

3.31 

3.41 

3.52 



Hi 

Fei 

Fen 

Fein 

Hn 

Om 

Cm 

Feiv 

Fev 

Orv 

Civ 

Niv 

Fevi 

Ov 

Nv 

Fe vn 

Cv 

Ovi 

Nvi 

Fe vm 

Arix 

Ovn 

Cvi 

Ne vm 

Arx 

Feix 

Nvn 

Fex 

Arxi 

Neix 

Fexi 

O vm 

Fexn 

Caxm 

Fexm 

Arxn 

Fexiv 

Fexv 

Sxn 

Fexvi 

Mgxi 

Fe xvn 

Caxiv 

Sxm 

Nex 

Fe xvm 

Sxrv 

Fexix 

Fexx 

Fexxi 

Fe xxn 

Fe xxm 

Fe xxiv 



< 1.00 
<0.63 

<0.03 

<0.25 

<0.03 

<0.25 

0.87 ± 0.23 



0.17 + 
<0. 

1.26 ± 
<0. 
<0. 

1.8 ± 
1.78 ± 
0.62 + 

<0. 
0.07 ± 

8.5 ± 

3.6 + 
0.3 + 

0.35 ± 
< 1. 

1.51 ± 

0.22 + 

<0. 

5.0 ± 

0.13 ± 

21.9: 

0.10 ± 
<0. 

0.26 ± 
<0. 

0.30 + 
<0. 

0.13 ± 

0.17 ± 

1.7 ± 
0.20 + 
0.13 ± 
0.24 + 

7.9 ± 
0.60 ± 
0.12 ± 
0.62 + 
0.48 ± 
0.28 + 

<0. 
<0. 
<2. 



0.20 
20 

0.15 

40 

06 

0.3 

0.22 

0.06 
25 

0.08 

0.6 

0.3 

0.3 

0.18 
26 

0.11 

0.04 
01 

0.6 

0.10 
= 2.1 

0.07 
05 

0.05 
00 

0.05 
05 

0.07 

0.17 

0.5 

0.04 

0.03 

0.07 

2.3 

0.07 

0.12 

0.06 

0.08 

0.10 
32 
40 
51 



1.20 
0.00 
0.00 
0.00 
604980 
0.09 
0.10 
0.00 
0.01 
0.86 
0.25 
0.09 
0.03 
1.24 
0.11 
0.09 
2.01 
1.76 
0.62 
0.17 
0.01 
8.55 
3.53 
0.38 
0.01 
0.09 
1.51 
0.20 
0.04 
5.00 
0.29 
21.97 
0.26 
0.18 
0.19 
0.06 
0.11 
0.05 
0.20 
0.04 
1.74 
0.27 
0.13 
0.09 
7.71 
0.45 
0.10 
0.66 
0.41 
0.30 
0.60 
1.39 
2.83 



0.1 
0.4 
0.1 



-0.7 
0.1 
0.0 

0.8 
-0.1 

0.4 
-0.4 

1.9 

0.0 
0.5 

0.0 

-1.6 
0.0 

-2.4 

1.7 

3.7 

-1.1 
0.8 
0.0 

-1.5 
0.0 
2.2 
0.1 
2.0 
0.2 
-0.7 
0.8 
-0.3 



79 




100 



99 
97 
100 
99 
90 
80 
85 
86 
41 
39 
28 
19 

7 

3 

3 





































































20 




1 
1 

3 


1 

10 
20 
15 
14 
59 
61 
72 
78 
91 
76 
94 
97 
38 
28 
34 
34 
32 
5 
2 
1 
1 


1 

























1 




8 











3 
2 
21 
3 
3 

60 
63 
66 
66 
68 
80 
98 
99 
92 
100 
75 
100 
96 
100 
98 
100 
99 
95 
95 
69 
59 
71 
61 
34 
11 
10 
1 













9 

















1 

6 






10 





7 


16 



4 



2 



1 

5 

5 

29 
41 
29 
39 
38 
89 
86 
98 
90 
35 

3 














81 



































3 







5 









8 

















2 






28 



3 

1 

10 
65 
97 
100 
100 



<0.08 
<0.50 

<0.32 

<0.16 

<0.05 

<0.08 

0.72 ±0.19 

0.1 ±0.3 

<0.16 

0.47 ± 0.24 

<0.01 

<0.03 

2.04 ± 0.25 

1.00 ±0.23 

0.59 ±0.10 

<0.20 

<0.03 

8.1 ± 1.0 

3.1 ±0.4 
0.22 ±0.12 

<0.50 
<0.04 
0.93 ±0.11 
<0.03 
<0.02 

2.2 ± 0.7 
0.05 ± 0.04 

13.5 ±2.4 
<0.03 
<0.05 

0.28 ±0.16 
<0.03 
<0.06 
<0.50 

0.17 ±0.07 

0.10 ±0.04 
<0.63 
<0.01 
<0.04 

0.09 ±0.14 
<0.63 
<0.08 
<0.13 
<0.50 

0.28 ±0.12 

<0.06 

0.2 ± 0.3 

<0.40 

<0.32 



<0.08 
<0.32 
<0.13 

<0.25 

<0.16 
<0.10 
<0.03 

<0.13 

<0.10 

1.0 ±0.5 

<0.50 

0.19 ±0.04 

<0.02 

0.7 ± 1.8 

0.06 ±0.10 

<0.04 

<0.10 

1.2 ±0.3 

0.62 ±0.11 
<0.03 
<0.05 
<0.03 

0.55 ±0.10 

0.26 ± 0.04 
<0.01 

2.0 ± 0.7 
<0.10 

5.1 ±1.5 
0.23 ± 0.05 

<0.02 

<0.06 

<0.01 

<0.05 

0.2 ± 0.4 

<0.04 

<0.03 

< 12.6 

0.19 ±0.04 

0.13 ±0.03 

0.16 ±0.09 

9.5 ± 1.4 

0.43 ± 0.07 

0.10 ±0.17 

0.71 ±0.05 

0.4 ± 0.3 

<0.06 

<0.16 

<0.40 

<0.79 



0.32 ± 0.25 



1.78 ± 1.08 



<0.32 



0.39 ±0.13 



270 ± 170 



160 ± 60 



- 10 ±50 
120 ± 60 

- 10 ± 50 



80 ± 20 
50 ± 20 



- 120 ± 40 

- 350 ± 220 

- 160 ± 40 

- 100 ± 20 

- 20 ± 190 

- 370 ± 400 

10 ± 110 

- 60 ± 230 

■ 640 ± 210 

■ 400 ± 130 

■ 400 ± 150 

- 230 ± 180 

- 270 ± 50 

■ 360 ± 1 10 

■210± 100 

■ 680 ± 300 

- 800 ± 220 



" Ionization parameter where the ion has its peak concentration in 10 v W m. 
* Element and ionization degree. 



'" Observed column density in 10"' 
d Predicted column density in 10 20 ' 



for model 1 . 
' m~ 2 for model 1 . 
e Difference in Ax between predicted and observed column density. 
' Percentage of ionic column density produced by component A (see Table[6}. 
s Percentage of ionic column density produced by component B (see Table[6]l. 
'' Percentage of ionic column density produced by component C (see Table[6]l. 
' Percentage of ionic column density produced by component D (see Table[6). 
' Percentage of ionic column density produced by component E (see Table[6]l. 

Observed column density in 10 20 irT 2 for model 2, velocity component 1 (v= — 13 ± 11 km s _1 ). The velocity broadening a = 125 ± 8 km s _1 . 
' Observed column density in 10 20 irr 2 for model 2, velocity component 2 (v = - 319 ± 14 km s" 1 ). The velocity broadening a = 107 ± 9 km s" 1 . 
'" Observed column density in 10 20 irT 2 for model 2, velocity component 3 (v = - 770 ± 109 km s ). The velocity broadening a = 160 ± 120 km s . 
" Outflow velocity for the individual ion in km s . Only ions with a solid measure of the column density are included (i.e. no upper limits). 
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3.4. Absorption measure distribution (AMD) 

There has been a debate in the literature about whether the ab- 
sorption measure distribution defined here as A(£) = dNu/d^ 
is a smooth distribution spanning several decades in t; (see 
e.g. ISteenbrugge et al.l 120051) or consi sts of a limited num- 
ber of discrete com ponents (see e.g. ICostantini et al.l l2007t 
iKrongold et aTEool . 

We tested both alternatives as follows. We considered the 
total ionic column densities derived from Model 1 (see Table 0J, 
regardless of their velocity, and fit them to a model using discrete 
components and to a model with a continuous distribution. We 
simultaneously solved for the abundances of the elements. 

For Model 1, there is some arbitrariness in the assignment 
of ions to the two velocity components, in particular near the 
ionization parameter log £ = 2.1 at the division, as for those ions 
both velocity components will contribute. Thus the observed col- 
umn densities near that division are higher than the column den- 
sity of a single velocity component. This could result in the in- 
troduction of spurious ionization components or artificially en- 
hanced abundances. We tried to also do this analysis with the 
results for each velocity component separately, but that is prob- 
lematic. For Model 2, the error bars on the column densities for 
individual velocity components are relatively high, because RGS 
only partially resolves the lines of each component. This then 
gives too much uncertainty to deduce conclusive results. 

The first model we tested is a discrete model: 



Table 6. Parameters for the discrete distribution. 



A(£) = y>^-6), 



(2) 



Where //, are the total hydrogen column densities of the m com- 
ponents with ionization parameter £,. 

From our runs with Cloudy we obtain for each ion j curves 
for the ion con centration f,(f) relative to hydrogen as a function 
of £, assuming lLodders & Palme (2009) abundances. Given a set 
of values for £, and //,, it is then straightforward to predict the 
ionic column densities Nf. 



DO 



(3) 



with Bj the abundance in solar units of the parent element of ion 
j. We solve this system by searching grids of models for differ- 
ent values of £,, and determine the corresponding best-fit column 
densities //, from a least-squares fit to the data. The abundances 
are solved iteratively. We start with solar abundances and solve 
for Hi, Then for each element we determine its best-fit abun- 
dance from a least squares fit of its ionic column densities to the 
predicted model of the last step. This procedure is repeated a few 
times and converges rapidly. 

It should be noted that since we do not measure hydro- 
gen lines, the hydrogen column densities that we derive are 
nominal values based on the assumption of on average solar 
metal-to-hydrogen abundance for the ions involved. In fact, 
we derive only accurate relative metal abundances. Truly abso- 
lute abundances should be derived using UV data, but we de- 
fer the discussion on ab undances to later papers of this series 
(Steenbrugge e t al.ll201 lL and Arav et al., in prep). 

It appears that we obtain the best solution if we take five 
ionization components into account. Adding a sixth component 
does not improve the fit significantly, and by deleting one, two 
or three components^ 2 increases by 6, 16 and 200, respectively 
(refitting in each case). Our best fit then has x 1 - 42.5 for 29 



Component log £° 



N' 



A 
B 
C 
D 
E 



-0.33 ± 0.49 0.23 ± 0.09 

0.71 ±0.12 0.84 ±0.10 
2.01 ± 0.02 4.8 ± 0.4 
2.79 ± 0.06 5.7 ± 0.9 
3.60 ± 0.27 54 ± 73 



" Ionization parameter in 10 5 W m. 
b Column density in units of 10 24 irT 



ions included in our fit. The predicted model is shown in Table 
|U together with the individual contributions A\j to x 1 f° r eacn 
ion (i.e., x 2 — 2 AV/)- Negative values for A%j indicate lower 
observed ionic column densities than the model and positive val- 
ues higher observed ionic columns. The best-fit parameters are 
shown in Table [6] 

We did not include upper limits in our fit, and we also ex- 
cluded the argon lines because the predicted model is well below 
the marginal "detections" of Arix and Arx. For further discus- 
sion, we also include predicted column densities for hydrogen 
and C in and C iv, although we cannot measure lines from these 
ions in the RGS band. 

Next we consider a continuous AMD. It is impossible to 
make no a priori assumptions for the shape of A(£), but we 
minimize this as follows. We assume that logA(£) is described 
by a cubic spline for log£ between -3 and 4 with grid points 
separated by 0.2 in log£. The use of logarithms guarantees 
that A(£) is non-negative; the spacing of 0.2 corresponds to the 
typical scale on which ion concentrations change (making it 
much smaller causes oversampling with unstable, oscillatory so- 
lutions), and the range in £ covers the ions that are detected in 
the spectrum. Free parameters of the model are the hydrogen col- 
umn densities H, at the grid points and the abundances. We sol ve 
for this system using a genetic algorithm (ICharbonneaul 1 9951) . 

We made 200 runs with the algorithm, and kept the 117 runs 
that resulted in v 2 < v 2 . +1 with y 2 . - 39.6 the best solu- 

^ ^ min ^ mm 

tion. In Fig. [9] we show the median of all these 1 17 good solu- 
tions. The figure shows two strong, isolated peaks at log£ = 2.0 
and 2.8, corresponding to components C and D of Table|6] At a 
higher ionization parameter (log £ > 3), the range of component 
E of Table [6] there is also some AMD, but the detailed struc- 
ture is essentially unknown: there is a wide spread between the 
individual solutions that are acceptable. At a lower ionization 
parameter (log£ < 1), there is also some AMD but again not a 
well-determined structure. A hint for the presence of component 
B is that the median of the acceptable solutions is closer to the 
upper limit in the range of £ between 0.4 - 0.8. 

As a final test, we extended the model with discrete com- 
ponents and searched how broad the discrete components are. 
Replacing the ^-function by a Gaussian in (0, we get an upper 
limit to the cr of the Gaussians of 0.06 and 0.13 in log£ for the 
components C and D, corresponding to a FWHM of 35 and 80%. 
For the other components, there is no useful constraint. 

4. Discussion 

4.1. Foreground and emission features 

This paper focuses on the properties of the ionized outflow in 
Mrk 509, but given the quality of the data, a full description and 
discussion of all the features present in the spectrum is beyond 
the scope of this paper and will be given in a series of subsequent 
publications. The only clear detections of narrow emission lines 
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Fig. 4. The data and the best fit to the RGS spectrum (model 2). The wavelengths shown are the observed wavelengths. 
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Fig. 5. RGS spectrum continued. 
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Fig. 7. RGS spectrum continued. 
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Fig. 8. RGS spectrum continued. 
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Fig. 9. Continuous absorption measure distribution plotted as 
%A(g) for Mrk 509. Shown is the median of the 117 runs with 
X 2 < X 1 ■ + 1 ■ Error bars correspond to the minimum and maxi- 
mum value of the AMD at each value of £ for those runs with ac- 
ceptable^ 2 . Note the two strong and isolated peaks at log £ = 2.0 
and 2.8, respectively. The inset shows our results on a logarith- 
mic scale for a broader range of ionization parameter; dotted 
lines connect the minimum and maximum values of all accept- 
able solutions. 



are the Ovu forbidden emission line at 22.101 A, the Ovu in- 
tercombination line at 21.802 A, and a narrow Neix forbidden 
emission line at 13.70 A. Th e fluxes are consistent with earlier 
observations dDetmers et alj |2010). We also searched for RRC 
features of the most prominent ions (C, N, and O) and found 
several possible weak RRCs. We do not detect any significant 
absorption due to neutral gas in the host galaxy. This means that 
we are observing the nucleus directly, which could indicate that 
we o bserve Mrk 509 almost face-on, as s uggested in previous pa - 
pers dPhillips et al.ll983hlKriss et alj2000tlKraemer et al.l2003l) . 

4.2. Outflow models 

The warm absorber in Mrk 509 (O vm column density of =* 2 x 
10 21 m~ 2 ) is deeper tha n the one seen in Mrk 279 (3 x 10 20 irT 2 , 
ICostantini et al.l 120071) but shallower with r espect to those in 
NGC 5548 for instance (3 xlO 22 nT 2 , see e.g. Steenbrugge et al. 
l2005h or NGC 3783 (4 x 10 22 nr 2 , see e.g. Beharet al. 2003). 
Nevertheless the high quality of this dataset allows for a thor- 
ough investigation of the outflow properties. While the main goal 
of the campaign is to localize the outflow, which requires inves- 
tigating the ten individual observations, the integrated 600 ks 
spectrum is crucial for a full description of the properties (such 
as outflow velocity and ionization structure) of the outflow. To 
obtain the most accurate information about the true structure of 
the outflow we now compare the different models of the outflow. 
The first model (Model 1) is a very simple description, with 
only one velocity component for each ion. The velocity dis- 
persion we obtain for the component that includes the Ovn 
and Oviii ions is 158 km s , This is larger than what is ob- 
tained from the curve of growth analysis using a single veloc- 
ity co mponent for these ions (96 and 113 km s , respectively, 
Kaastra et al. 201 lb). The reason for this difference is that hav- 
ing only one velocity component for these lines is an oversim- 
plification. Adding a second velocity component for all ions 



(Model 2), improves the fit of the strong oxygen lines (O vn and 
O vm), with the sum of the velocity dispersions larger than for 
the single componen t case, but the total column density similar 
dKaastra et al.ll2.Ql Tbl) . Thus the total ionic column densities for 
models 1 and 2 are consistent with each other. 

To compare the slab models with model 3 we first need to 
convert the ionic column densities we measure into an equivalent 
hydrogen column density. There are two ways to do this. One 
is to assume that every ion occurs at the ionization parameter 
where its concentration peaks as a function of log ^. This holds 
for some ions, but for others there is a wide range of ionization 
parameters where the ion makes a significant contribution. The 
alternative is that we take the full AMD method described in 
Sect. 13.41 It is useful to compare these two methods so that we 
can see if there are major differences in the results and if these 
possible differences affect our conclusions. The results for the 
first method using Model 2 are shown in Fig. [TUJ for both the 
slow and fast velocity components. Only ions for which we have 
a significant column density measurement are shown. The results 
for the AMD method are shown in Fig. [9] What can be seen is 
that the AMD method clearly shows a discrete distribution of 
column density as a function of the ionization parameter. There 
is a clear minimum between the peaks at log £, - 2.0 and 2.8, 
where the column density is more than an order of magnitude 
less than at the two surrounding peaks. This indicates that there 
is almost no gas present at those intermediate ionization states. 
The simplified method (the one where we assume that every ion 
occurs at a single £, value) does show enhancements near the 
mean peaks of log £ = 2 and 2.8, but there are no clear minima 
in the distribution, although for the fast component there seems 
to be some hint for a minimum near log £ = 0.5. What is clear 
from this comparison is that the simplified method is unable to 
uncover essential details in the AMD. This is because not all ions 
are found at their peak ionization parameters. 

Another main difference between the models is that the 
slab models (in contrast to Model 3) yield completely model- 
independent ionic column densities (i.e. no SED or ionization 
balance or abundances are assumed). This is an advantage if the 
atomic data for certain ions are uncertain, as the fit will not try 
to correct for this by changing the overall fit parameters or by 
poorly fitting this particular ion. From the measured ionic col- 
umn densities we then can obtain the distribution of total hydro- 
gen column density as a function of the ionization parameters. 
However at this step it requires the input of an ionization bal- 
ance, hence an SED. 

Model 3 is a direct fit using the ionization balance to predict 
the ionic column densities. The advantage of this model com- 
pared to Models 1 and 2 is that all ions, including those with a 
small column density are taken into account, and all the ionic 
column densities are connected through a physical model. There 
are fewer free parameters, so in principle a more accurate de- 
termination of the following parameters: Nh, £, <T and v can be 
obtained. The AMD method is a good combination of both mod- 
els because the ionic column densities are determined model in- 
dependently by the slab models. Then the ionization balance is 
used to produce an AMD and obtain the number and parameters 
of xabs components that are needed to describe the data properly. 

Model 3 and the AMD method make use of the ionization 
balance as determined from the SED. The resulting stability 
curve for the photoionized gas is shown in Fig.QT] Components 
with the same S (in units of ) are in pressure equilibrium. Here 
S is defined in the following way: 



E = 9610x£/7\ 



(4) 
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Fig. 10. The derived hydrogen column density for every de- 
tected ion (see Table @). We added archival UV data for C in, 
C iv, N v, and O vi for comparison (shown in red). The top figure 
shows the distribution for the slow component, the bottom one 
shows the same for the fast component. 



where T is the temperature and £ is the ionization parameter in 
1CT 9 W m. The sections of the curve with a negative slope are un- 
stable to perturbations. Not all components appear to be in pres- 
sure balance. For both velocity components the low-ionization 
gas (A2, Bl) is not in pressure equilibrium with the higher ion- 
ization gas. This could indicate that the different gas phases are 
not colocated or that other forces (i.e. magnetic) are involved 
to maintain pressure equilibrium. It appears that most of the 
outflows in Seyfert 1 galaxies show gaps in the AMD (Behar 
2009 ). This is proba bly due to (thermal) instabilities in the gas 
dHolczer et al.ll2007l) . The exact nature of these apparent insta- 
bilities is still unclear (although a thermal scenario indeed seems 
plausible at the moment, based on the cooling curves, such as 
shown in Fig.fTTt. 



4.3. Structure of the outflow 

Much work has already been done investigati ng the structure of 
the ionized outflows in other AGN (see e.g. Steenbrugge et al.l 
I2005L iHolczer et al.ll2007t ICostantini et alJl2007l for some ex- 
amples). In most cases a wide range of ionization states has been 
detected, sketching the picture of a continuous distribution of the 
hydrogen column density as a function of £. However, there are 
also indications of a lack of ions in a certain temperature regime, 
where the ionized gas is in an unstable region of the cooling 




Fig. 11. The cooling curve derived for the SED with the ioniza- 
tion parameters obtained from Model 3 over-plotted. The circles 
indicate the ionization parameters for the slow velocity compo- 
nent, while the squares shown are for the fast outflow. Regions 
where heating (H) or cooling (C) dominate are also indicated. 



curve. In Mrk 509 such unstable regions occur for log f between 
2.4 and 2.8 and between 3.5 and 4. We have determined that the 
outflow in Mrk 509 is not continuous, but has discrete compo- 
nents, at least in the range of log £ = 2 - 3. There are two main 
components, one at log £ = 2.0 and one at 2.8. We also see a 
clear trend toward increasing column density for higher ioniza- 
tion states. 

We first discuss the ionization structure. The most pro- 
nounced component in our spectrum is component C (Table [6). 
The ionization parameter and total column density derived from 
our slab fit (Table [6), as well as the direct xabs fit (Table |5), are 
fully consistent. From the analysis in Sect. I3.4l we find that this 
component is discrete and spans a very narrow range in ioniza- 
tion parameter: the FWHM is 35%. Interestingly, according to 
our model (Table |4]i component C contributes 50% or more to 
the total ionic column density of 17 of the detected ions in our 
data set. Those ions span a range of log£ = 1.15 (O vn) to 2.42 
(S xm) in ionization parameter (see first column of Table©. 

The next most important component is component D at 
log£ = 2.79. It is responsible for the more highly ionized iron 
(up to Fexx) and the sulfur ions. It is mainly visible in the 
high-velocity component. Again the direct xabs fit and the de- 
rived ionization parameter are fully consistent with each other, 
but the column density obtained from the xabs fit is smaller by 
a factor of 3. This could be due to xabs component C2 (log £ = 
2.2), which also produces ions present in component D (log t; = 
2.79). The other important component is B, which is responsible 
for most of the lower ionized carbon, nitrogen and oxygen ions. 
The ionization parameter and column density are fully consistent 
for both the derived model and the direct xabs fit. 

Components A and E are also fully consistent with our 
Model 3, however they are only based on a few ions, so their 
exact column densities and ionization parameters are uncertain. 
Especially for the low-ionized gas, the lack of a strong UTA in 
Mrk 509 means that we only have upper limits on the column 
densities of the low ionized iron ions, up to Fe x. This is also 
why the AMD distribution shown in Fig. |9]has large uncertain- 
ties below log £, — 0. 
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Earlier observations have detected only a few of the five 
components shown in Table[6] owing to the poorer quality of the 
data. Using the Chandra HETGS, Yaqoob et al. (2003) detected 
mainly component C, since the sensitivity of the HETGS is lim- 
ited at long er wavelengths a nd component C is the strongest 
component. ISmith et alj d2007l) detected mainly component B2, 
C (mixture CI and C2), and D2 (possibly blended with E2). The 
outflow velocities in their analysis are different, as mentioned 
before, and they observe an inverse correlation between the out- 
flow velocity and the ionization paramete r. This analysis was 
based on the 2000 and 2001 archival data. iPetmers et all d2010h 
analyzed earlier archival data (2005 and 2006) of Mrk 509 and 
found component B (possible mixture of the velocity compo- 
nents), CI, and D2 (possibly blended with E2). It is clear from 
these comparisons that, although all these earlier observations 
detected the main components of the outflow, in order to obtain 
a more complete picture of the outflow, we need a high-quality 
spectrum, like the one shown in this paper. 

Two of the three velocity components that we d etected are 
consistent with earlier results, including the UV data (iKriss et al.l 
120001: lYaqoob et al.ll2003t iKraemer et ai]|2003l) . and they corre- 
spond to the two main groups of UV velocity components, one 
at systemic velocity and the other at - 370 km s . Also in the 
X-ray regime, there is evidence of multiple ionization states for 
the same outflow velocity, such as components C2 and D2 in 
Table [5] Also there are components that show a similar ioniza- 
tion state, but different outflow velocity, i.e. components CI and 
C2 in Table[5] Due to the almost zero outflow velocity of compo- 
nents B 1 and C 1 , one could argue that these may be related to the 
ISM of the host galaxy. The UV spectra, with their much higher 
spectral resolution, can unravel the outflow, ISM, and redshifted 
high-velocity clouds (Kriss et al. 201 1). Generally speaking, the 
ionization parameter of the UV components is much lower than 
those of the X-ray components detected here. This could indicate 
that the UV and X-ray absorbers are cospatial, but have different 
densities. A full discussion of the connection between the UV 
and X-ray absorbers , as well as the geo metry of the absorber, 
will be presented in lEbrero et al.l (1201 lh . where the simultane- 
ous HST COS and Chandra LETGS data will be compared. 
We do not clearly detect the 200 km s redshifted compo- 
ne nt, which w as found in the UV data (velocity component 7 
of IKriss e t al. 2000). There is some indication that there could 
be an O vi component at that velocity. There is some extra ab- 
sorption at the red side of the line in Fig [6] at 22.78 A. We only 
obtain an upper limit for the O vi column density in this veloc- 
ity component of 10 20 irT 2 . However, this is consistent with the 
lower limit from the UV data, which is 10 19 irT 2 . 

The highest velocity outflow component is only significantly 
detected in two ions (Mg xi and Fe xxi). FigureQ~2]shows the two 
absorption lines fitted with a - 770 km s -1 velocity and a - 300 
km s _1 velocity (just as component two in model 2). The im- 
provement using the - 770 km s component is Ay 2 = 16. We 
checked whether this component is also detected in other ions, 
but most of them (apart from Ne ix) only yield upper limits to 
the ionic column density. We checked that the line profiles are 
the same for the separate RGS 1 and RGS 2 spectra and also 
for the first and second order spectra. In all cases the line pro- 
file is consistent with a 770 km s blueshift. Also a possible 
Mg xi forbidden emission line cannot play a role here, because 
it is too far away (about 0.3 A) to influence the line profile in 
any way. This velocity component is consistent with an earlier 
Chandra HETGS observation , where there w as an indication of 
this velocity component dYaqoob et al]l2003l) . A proper expla- 
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Fig. 12. The absorption lines of Mg xi and Fe xxi, which show a 
blueshift of 770 km s -1 . The dashed-dotted line shows a model 
with a blueshift of 300 km s and the solid line shows the model 
with a blueshift of 770 km s" 1 . 



nation of why this component is only clearly detected in these 
two ions and not in other ions with a similar ionization param- 
eter is currently lacking. A trend visible in Tableland Fig. Qj] 
is that the higher ionized ions have a higher outflow velocity. 
Fitting a constant outflow velocity to the data yields a value of 
70 + 9 km s with^ 2 = 72 for 24 d.o.f. A linear fit to the points 
gives a slope of 0.62 + 0.07 and improves the y 1 to 46 for 24 
d.o.f. If we instead fit a powerlaw, a relation of v =* ^°- 64±01 ° is 
obtai ned, with a to t al y 2 of 34 for 24 d.o.f. The MHD models 
of Fukumuraet al. (2010) predict v =* £° 5 , which is consistent 
with the relation found here. However, it has to be noted that 
due to the blending of multiple velocity components (including 
gas that might not be outflowing at all), additional uncertainties 
are introduced that could affect the results, selecting only those 
ions, which clearly show blueshifts as well as more accurate out- 
flow velocities, would be needed to investigate this trend further. 
The O iv ion shows a large redshifted velocity, undetected in the 
other ions. Most likely this is due to the blending of the Oiv 
absorption line with the O i line from the Galactic ISM at z = 0. 
This blend ma kes it difficult to det ermine the centroid of the O iv 
line exactly dKaastra et al.ll2011bh . 
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Fig. 13. Outflow velocity vs. £. A general trend of increasing 
outflow velocity for ions with a higher ionization parameter can 
be seen. A negative velocity indicates outflow. 



4.4. Density profile 

A recent study has used the observed AMD to construct t he ra- 
dial d ensity profile of the outflow in a number of sources (Behar 
2009). Such an analysis is justified as long as the AMD is a 
smooth, continuous function of £, with the possible exception of 
unstable branches of the cooling curve, where gas may disappear 
to cooler or hotter phases. However in our case such an analy- 
sis is not justified, at least not for the range of log £, between 
~ 2 and 3 (components C and D). Here we clearly see narrow 
peaks in the AMD. At least for these components, this hints at 
rather localized regions with a limited density range, rather than 
to a large-scale outflow. At lower ionization parameters (com- 
ponents A and B), we cannot exclude a continuous distribution, 
owing to the limitations imposed by the line detection from the 
relevant ions. For these components, the situation is more com- 
plex due to the presence of both higher and lower velocity gas. 
Similarly, based on our analysis we cannot distinguish whether 
component E has a single component or a broader distribution 
on the second stable branch of the cooling curve (Fig. fTTT i. Gas 
on the third stable branch, at a very high ionization parameter, 
escapes our detection completely because of the lack of suit- 
able lines in the RGS band. At best it could show strong lines 
from Fexxv or Fexxvi in the Fe-K band near 6.7 - 7.0 keV, 
but the limited spectral resolution of EPIC combined with the 
likely moderate column densities prohibit us from detecting such 
a component in our data. 



absorption measure distribution, but instead shows two strong, 
discrete peaks. At the highest and lowest ionization parameters, 
we cannot distinguish between smooth and discrete components. 
We also have found indications of an increasing outflow velocity 
versus ionization parameter. Large, dedicated multiwavelength 
campaigns such as this are the way forward, as this is currently 
the best method to investigate and characterize the outflows in 
the local Seyfert galaxies. 
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5. Conclusions 

We have presented one of the highest signal-to-noise RGS spec- 
tra of an AGN. With the almost unprecedented detail in this 
dataset, we could detect multiple absorption systems. The ion- 
ized absorber of Mrk 509 shows three velocity components, 
one at -13 ±11 km s" 1 , one at -319 km s , and a tentative 
high-velocity component at - 770 km s . The first two compo- 
nents are consistent with the main absorption troughs in the UV. 
Thanks to the high-quality spectrum and the accurate column 
densities obtained for all ions, for the first time it has been shown 
clearly that the outflow in Mrk 509 in the important range of 
log £, between 1-3 cannot be described by a smooth, continuous 
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Table A.l. Updated line list for important ions in the Mrk 509 
spectrum. 



Ion 


Wavelength (new) Wavelength (old) 


Reference 1 




(A) 


(A) 




Nvi 


28.7875 


28.7870 


1 


Ora 


22.9400 


22.9784 


2 


Ora 


23.0280 


23.0489 


2 


Ora 


23.0710 


23.1092 


2 


Ov 


19.3570 


19.3251 


2 


Ov 


19.9680 


19.9242 


2 


Ovi 


18.2699 


18.2896 


2 


Ovi 


18.2700 


18.2897 


2 


Ovi 


18.5869 


18.6059 


2 


Ovi 


18.5870 


18.6060 


2 


Ovi 


19.1798 


19.1355 


2 


Ovi 


19.1805 


19.1362 


2 


Ovi 


19.3789 


19.3412 


2 


Ovi 


19.3791 


19.3414 


2 


Ovi 


22.0189 


22.0063 


3 


Ovi 


22.0205 


22.0079 


3 


Ovii 


17.7683 


17.7680 


1 


Ovu 


18.6284 


18.6288 


1 


Ne vm 


13.6533 


13.6460 


4 


Nevm 


13.6553 


13.6480 


4 


Sxm 


32.2380 


32.2420 


5 


Sxiv 


30.4330 


30.4270 


5 


Sxiv 


30.4750 


30.4690 


5 


Fe xvn 


15.2610 


15.2650 


6 


Fe xvm 


14.3720 


14.3780 


7 


Fe xvm 


14.5340 


14.5400 


7 


Fe xvm 


14.5710 


14.5550 


7 


Fexix 


13.4620 


13.4650 


7 


Fexix 


13.5180 


13.5210 


7 


Fexix 


13.7950 


13.7980 


7 


Fexx 


12.8240 


12.8130 


7 


Fexx 


12.8460 


12.8270 


7 


Fexx 


12.8640 


12.8470 


7 


Fexx 


12.9150 


12.9040 


7 


Fexxi 


12.2840 


12.2860 


7 
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(7) Brown et al. (2002). 
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Appendix A: Improving the atomic data 

Due to the high quality of the Mrk 509 dataset, we have also 
updated some of the laboratory wavelengths of the important 
ions detected in the Mrk 509 spectrum. Table IA.l1 shows all the 
updated lines. The Oin lines are n ot resolved, so the strongest 
line at 23.071 was used CGui 2006) and the other two lines were 
shifted by the same amount. 



